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Shock waves are supersonic disturbances propagating in a fluid and giving rise to dissipation 
and drag. Weak shocks, i.e., those of small amplitude, can be well described within the 
hydrodynamic approximation. On the other hand, strong shocks are discontinuous within 
hydrodynamics and therefore probe the microscopies of the theory. In this paper we consider 
the case of the strongly coupled M = 4 plasma whose microscopic description, applicable for 
scales smaller than the inverse temperature, is given in terms of gravity in an asymptotically 
AdS 5 space. In the gravity approximation, weak and strong shocks should be described 
by smooth metrics with no discontinuities. For weak shocks we find the dual metric in a 
derivative expansion and for strong shocks we use linearized gravity to find the exponential 
tail that determines the width of the shock. In particular we find that, when the velocity 
of the fluid relative to the shock approaches the speed of light v — > 1 the penetration depth 
£ scales as £ ~ (1 — v 2 ) 1 ^ 4 . We compare the results with second order hydrodynamics and 
the Israel- Stewart approximation. Although they all agree in the hydrodynamic regime of 
weak shocks, we show that there is not even qualitative agreement for strong shocks. For 
the gravity side, the existence of shock waves implies that there are disturbances of constant 
shape propagating on the horizon of the dual black holes. 
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1 Introduction and summary 



The AdS/CFT correspondence [1, 2, 3, 4] realizes in practice the idea of describing gauge 
theories in terms of strings [5] and provides a new tool to study gauge theories at strong 
coupling. Recently, motivated by experimental and theoretical work on heavy ion collisions 
at RHIC, there was a wave of interest in using AdS/CFT to study the dynamics of strongly 
coupled plasmas, for a review see [6, 7, 8, 9, 10] and references therein. The main initial 
focus was on the low viscosity of the fluid and the drag force experienced by a quark moving 
through the plasma. More generically, in conformal plasmas, such as the one studied in 
AdS/CFT, the hydrodynamic description of the theory is valid up to length scales of order 
of the inverse temperature. Below that scale standard hydrodynamics breaks down, and we 
need to resort to the dual gravity description, which does not break down until the (much 
smaller) bulk string scale is reached. Alternative descriptions that are not based directly 
on microscopies, such as the Israel-Stewart theory, can be tested by comparison with the 
microscopic theory provided by the dual gravity. 

In the present paper, we focus on a particular phenomenon in fluid dynamics, namely, 
shock waves. These are supersonic disturbances that propagate in the fluid and are typically 
produced by an object moving supersonically in the fluid or by a localized release of energy 
such as in explosions or collisions. In ideal fluids they are described by a surface of discon- 
tinuity where the normal velocity and the pressure have a jump. In the frame where the 
shock is at rest, the fluid goes from supersonic to subsonic, with the kinetic energy of the 
fluid converted into pressure and heat. This process is irreversible and generates entropy. 
In fact, it is the only dissipative mechanism in the case of strictly ideal fluids. In practice, 
of course, the fluid is never ideal so at small scales the viscosity also causes dissipation. If 
the shock wave is produced by an object moving in the fluid it also generates drag. Away 
from the shock the hydrodynamic approximation is good and allows to compute the jumps 
in velocity and pressure by using the conservation of energy (and other charges) across the 
shock. However in the region of discontinuity the gradients are in general large and a mi- 
croscopic theory is necessary to have an appropriate description. In fact this breakdown of 
hydrodynamics is the reason why dissipation occurs in shock waves even for ideal fluids. For 
that reason, shock waves are not only phenomenologically interesting but, on the theoretical 
side, can be considered as a useful probe of the system. In practice, they are usually studied 
in non-relativistic systems and in fluids such as air and water. The relativistic case is of 
interest in astrophysics and in particular in heavy ion physics, which is the system closest to 
the one studied in this paper. 
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In the context of heavy ion physics there are different situations where shock waves can 
appear. One is during the period of the creation of the Quark Gluon Plasma. A simple model 
of particle production which gives reasonable results is Landau's hydrodynamic model [11]. 
This model assumes that the fluid thermalizes instantaneously and then evolves according 
to ideal hydrodynamics. All the entropy is created at the initial stage as a result of the 
thermalization process. In a purely ideal hydrodynamical model the entropy can be generated 
only from shock waves which would therefore play an important role. However, initially the 
quark gluon plasma is far from equilibrium and therefore a hydrodynamical approximation 
might not be valid. 

Within AdS/CFT entropy production in a heavy ion collision has been studied in [12] 
where two shock waves moving at the speed of light collide creating trapped surfaces. Also, 
[13] study the problem of reconstructing the bulk metric from the boundary. Several other 
aspects are examined in [14, 15]. Related work on gravitational shock waves appears in 
[16, 17, 18, 19, 20, 21, 22, 23, 24, 25]. 

Another instance where shock waves might appear in heavy ion collisions is the process 
where a hard parton moves through the plasma. The hard parton can be a heavy quark, a 
meson or a gluon. Calculations in AdS/CFT [26, 27, 28, 29, 30] show that there is a Mach 
cone created. That is, the energy density of the wave is concentrated in a angle close to 
the Mach angle. This is highly suggestive of the formation of a shock wave that cannot be 
described by the linearized approximations used in previous calculations. The real space 
profile of the stress energy tensor of a moving quark, calculated in [27] also supports the 
idea. Independently from AdS/CFT the possibility of shock waves in heavy ion collisions 
has been explored both theoretically and experimentally in [31, 32, 33, 34, 35, 36, 37, 38, 39, 
40, 41, 42, 43, 44, 45, 46, 47, 48]. 

2 Shock waves in hydrodynamics 

Hydrodynamics can be used to study the region away form the shock where the gradients 
are small. For weak shocks, i.e., those of small amplitude, a hydrodynamic approximation 
that includes viscosity resolves the discontinuity and provides a smooth description of the 
shock wave. This can be corrected at higher orders in gradients if so desired. In this section 
we start by considering the ideal fluid case for arbitrary shocks and then the higher order 
hydrodynamic approximation for weak shocks. We also analyze the Israel- Stewart theory to 
compare later with the microscopic results. 
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Figure 1: A supersonic flow reaches a fixed object creating a stationary shock wave. In heavy 
ion physics this scenario might be realized by a heavy parton moving through the plasma as 
suggested by the existence of a Mach cone. 

2.1 Ideal Hydrodynamics 

Ideal hydrodynamics is valid far from the shock and determines the main properties of shock 
waves. These are well known [49] and we describe them in this subsection as applied to our 
particular fluid. In relativistic ideal hydrodynamics the energy momentum tensor is given 
by 

T» v =prf w +(p + p)ui t u v , (1) 

where p is the pressure, p the energy density and w M the four velocity of the fluid. The 
equations of motion are given by conservation of T M ", 

d^ v = , (2) 

and the equation of state, which for the strongly coupled M = 4 plasma is dictated by 
conformality and reads 1 

p = 3p = 3(ttT) 4 . (3) 
J We follow the conventions of [50] with regard to the normalization of the stress energy tensor. That is 

o 2 

our stress tensor is related to the conventionally defined one by Tg£ r = j^-T^ . In gravity this is reflected 

[51 [51 

in TJ£L = 16wG JT^ where G J is the five dimensional Newton's constant. 
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For a shock wave with a planar front, located for convenience at x = 0, both T tx and T xx 
are constant throughout the fluid: 



T tx = ApuV = 4p- , (4) 

1 — v 2 

T xx = p (1 + 4u") = • (5) 

1 — v z 

The possibility of a shock wave arises because there are different values of pressure and 
velocity that give the same value of T tx and T xx . For example, the ratio T tx /T xx depends 
only on the velocity and, as is seen from the plot in fig. (3) different velocities can give the 
same ratio T tx /T xx . More precisely, if we take that for x — > — oo, v — V\ and T = T\ and 
for x — > oo, v = v 2 and T = T 2 we can ensure that T tx and T xx are constant by imposing 

1 9v 2 - 1 „ „ f 9vj -P 1/4 



V2 = ^— , P2= PiiTTT 2Y ' T 2 = T 1 [ —± 2V • 

3vi 3(1 - vf) \3(1 - vf) J 

A particular case is V\ — v 2 = v s where v s = 4^ is the speed of sound. In that case the 
jumps in v and T vanish. In general we have 

1> B1 >-L>„ 2 >I (7) 

as can be seen from the relation V\V 2 = |- In particular, if v 2 = 1/3 on one side of the 
shock, the fluid on the other side moves at the speed of light. In principle, the conservation 
laws allow one to switch the values of the velocity between the front and back of the shock, 
i.e., take V\ < v s and v 2 > v s , but such solution would convert thermal energy into kinetic 
energy violating the second law of thermodynamics. In figure (4) we show the ratio of the 
pressures on the two sides. The supersonic side of the shock has a lower pressure which 
goes to zero when v\ approaches one. For an object moving through the fluid this has the 
effect of changing the pressure that the object experiences. Indeed, defining the projector 
perpendicular to u M as P^ u = jf" + u^u u we obtain, from eq.(l): 

P^d a T av = 0, d lt T + u v d v (Tu li ) = (8) 

For a stationary solution, taking the \i = component we obtain: 

u% ( 7 T) = (9) 
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Figure 2: Sketch of a shock wave in the rest frame of the interface. For ideal hydrodynamics 
there is a discontinuity at x — 0. Including the higher order terms in the expression of the 
stress energy tensor resolves the discontinuity The conventions of the paper are that the 
fluid moves to the right. The left hand side is supersonic and the right hand side is subsonic 
since the opposite configuration violates the second law of thermodynamics. 

after replacing u° = 7, where 7 is the Lorentz factor. We see that 7T is conserved along a 
streamline. In particular for a situation as in fig. 1 we can follow a streamline from infinity 
to the point A and obtain, with and without a shock wave, the pressure: 

(Shock) 4 27l£ p >° Sh ° ck ) _ „4 _ Poo {m 

Pa -I2P2- ( 9 ^_i)(i_^)^' Pa -TooPoo- {1 _ v 2 q) 2 > ^ 

where we have used eq.(3) and the matching conditions (6) for the pressure and velocity 
across the shock wave. 

Another interesting aspect of shock waves is the entropy production that is associated 
with them. In ideal hydrodynamics the entropy current is given by 



s 



-4tt 4 TV, (li; 



where we have chosen the normalization to give the entropy density of a fluid at rest. The 
difference between s° on the two sides of the shock is given by 



As — s° subsonic s° supersonic — Air ( T2\/l + «2 1 + u\ 



" vr^v Hi--? 




(12) 
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Figure 3: The ratio T tx /T xx for different velocities. For V\ > | the equation T tx /T xx = const. 
has two solutions and hence allows for a jump in velocity. When v\ — ^= there is no jump. 

The entropy production comes from the difference of the two fluxes 
A** = -s x = 4?r27? (s-^ ( 9V * ~ 1 

° subsonic ° supersonic r 9 I \ -1 2 

V 1 - "1 \ V 1 - w i 

Kinetic energy from the supersonic side is transformed into thermal energy in the subsonic 
side and this creates entropy. Notice also that the ideal hydro is correct far from the shock, 
so this calculation gives the correct entropy production even if the fluid is not ideal. 

2.2 Viscous Hydrodynamics 

The stress-energy tensor for relativistic hydrodynamics can be organized in a series expansion 
in powers of jjf where L is a typical length scale over which the four-velocity changes and 
T is the temperature. In such an expansion the first order term is the viscous term. Since 
the plasma we are interested in is conformal, its bulk viscosity is zero. The shear viscosity 
is given by ^ = [51]. The stress energy tensor to the first order is given by 

= (vrT) 4 {rf v + Au"u v ) - 2(nT) 3 (r^ , (14) 

where 

pHV = ^ + u a u u ^ a uv = p^p^ d{aUp) _ ^P^d a U a . (15) 
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Figure 4: The ratio of pressure for the two sides of a shock wave as a function of the incident 
velocity V\. When v\ < 1/3 there is no shock. For V\ < ^,t ; 2 > the pressure and the 
temperature on the supersonic side (2) are smaller than those on the subsonic side. Kinetic 
energy from the supersonic side is transformed to thermal energy on the subsonic side thereby 
increasing the entropy. The entropy production in this process is discussed further in section 
(2.3). 

Notice that 

T^u u = -3(ttT)V . (16) 

which can be taken as the definition of u M and T, namely is a time like eigenvector of T^ u 
whose eigenvalue is — 3(?tT) 4 . This definition can be used at any order in the hydrodynamic 
expansion and is sometimes called the Landau frame. 

Now we would like to see how a weak shock wave is resolved if the effects of viscosity are 
included. We consider a flow where the four-velocity and temperature are functions only of 
x: 

T = T{x), u^ = (-v/1 + M (x) 2 ,m(x),0,0) . (17) 

As before, conservation of the energy-momentum tensor implies that the components T tx 
and T xx are constant throughout the fluid. They are now given by 

T tx = 4( 7 rT) 4 u(a;)v / l + u{x) 2 - \{>KTfu{x)^l + u(x) 2 u'(x) = f tx , (18) 

T xx = (ttT) 4 (1 + Au{xf) - ^(vrT) 3 (l + u{xf)u\x) = f xx . (19) 
The asymptotic behavior of T,u(x) determines the constants T tx ,T xx . Conversely, in terms 
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of T tx ,T xx and u(x) the temperature is given by 

*r=^(Wi±^-f-| , (20) 

which is in fact valid to all orders in the hydrodynamic expansion since it follows from the 
definition (16). 

Let us suppose that asymptotically on the supersonic side the temperature and velocity 
approach constants: 

Ti = lim T(x) = T( ) (1 - ^u^) , m = -^= + = lim u(x) . (21) 

X— > — oo \ <J / V 2 X— > — oo 

Four-velocity corresponds to the speed of sound in a conformal plasma. The remaining 
equation T tx = T tx gives a differential equation for u{x). Since first-order hydrodynamics is 
valid only for weak shocks we expand in a power series in Uoo. It then becomes clear that it 
is useful to define a new variable 

e = 4 ^"°°x , (22) 
3 

in terms of which we have 

u(0 = + (0 + «L^(2) (0 + • • • • (23) 

The equations of motion imply that <Sitm(£) satisfies the equation (primes denote derivation 
with respect to £) 

= 5^) - 1 , (24) 

with solution 

<J«(!) = tanh(-0 . (25) 

This has the same form as the solution for a weak shock in nonrelativistic hydrodynamics 
[49], and indeed Eq. (24) coincides with the first integral of the Burgers equation, familiar 
in that context. 

As we already noted, in ideal hydrodynamics one can freely exchange the two sides of 
the shock, so that in the rest frame of the shock the fluid's velocity may change either from 
subsonic to supersonic or from supersonic to subsonic. However, the existence of friction in 
the first order hydrodynamics breaks this symmetry and only the latter solution is allowed. 

The approach to the asymptotic values of T and u is described by T(x) ~ T asympt _ + 
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e iqx 5T, u{x) ~ u asymp t. + e tqx 5u, where T asympt is either 7\ or T 2 , depending on which of the 
asymptotic regions we are looking at. Expanding the equations for conservation of energy 
and momentum to the first order in ST, Su provides us with a system of two equations with 
two unknowns, 5T, 5u. Demanding that there is a non-zero solution determines q to be 2 



4 = ^HZ(3, 2 -1) (26) 

7Ti V 

where we use u asympt . = ^=5 ■ For weak shocks, we expand (26) around the speed of sound 
v ~ -j= + Sv to obtain 

4 = , (27) 

711 

which agrees with the explicit solution (25). For strong shocks, for which |g| is not small 
in comparison with T, there is no reason to expect (26) to be a good approximation. We 
compare it with other approximations in subsequent sections. 

2.3 Second order hydrodynamics and Israel- Stewart theory 

Let us now consider how shock waves are resolved in second order hydrodynamics and in 
Israel- Stewart theory. For the M = 4 plasma, the stress energy tensor has been computed 
to second order in [52, 50] (see also [53]), and is given by 



T* 4 " =(ttT) 4 {rf v + Au"u v ) - 2{nT) 3 a^ 

+ (ttT) 2 ((In 2)1% + 2T% + (2 - In 2) Ql^ + T% + T%) ) 



(2? 



where 



a*" = P» a P vP d {a u p) - l -P^ u d a u a , V = u a d a (29) 



TZ = e^a^ujp , Tg = d a uV , (30) 

T 2 7 = o^al - l -P^a^a aP , (31) 

= V»Vu v - -P^Vu a Vu a , (32) 
3 

Tg = P^P^V (d (a u p) ) - l -P^P^V (d a u ) , (33) 



1 Notice that this gives real exponential that decay away from the shock. 



10 



Ifi = e aMll u a dPv? , (34) 

We follow the conventions of [50] where €0123 = 1 and the brackets denote symmetrization. 
The T 0x and T xx components of the stress tensor are again constant but now given by 



T 0x =u(x)^/l + u(xy[ 4(vrT(x)) 4 - ^(ttT(x)) 3 + 
(ttT(x)) 2 ( 3(4 - ln2)u'(a;) a + \u{x)u"{x) 



T xx =(vrT(a;)) 4 (l + Au{xf) - Ul + u(x) 2 )(irT(x)) 3 + 

(ttT(x)) 2 (1 + u(xf) (1(4 - ln2)u\x) 2 + \u{x)u"{x) 



(35) 



(36) 



First, we carry out linear analysis near the asymptotics at x — > ±00, where we expect 

T[xj T aS y m pi + 5Te iqx , u(x) = u asympt . + 5ue iqx . (37) 
Keeping only linear terms in ST, 5u we solve 

rptx rptx rpxx rpxx (38) 

where T tx ,T XX are given by the asymptotic values 



T — 4:(7rT aS y mp t ) U asyrnp t. Y 1 + u asympt. i T ~ (^^asympt . ) (l + ^asympt .) ■ (^9) 

The 2 by 2 linear system for Su, ST has a nonzero solution only if its determinant is zero. 
This condition determines q to be 

%q = v / l^(l-V / l-2(3^-l)(2-log"(2))) (4Q) 
nT aS y mpt . v(2 - log(2)) 

where v = , n ° Si/ 2 mpt, = . This is intended as an improvement on the first-order formula (26). 
Note that the argument of the square root becomes negative for velocities greater than 
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This indicates that, not unexpectedly, computations in second order hydrodynamics should 
not be trusted beyond the weak shock regime, i.e., beyond velocities close to the speed of 
sound, \v — v s \ <C v s . Note that the velocity (41) is different from the velocity of discontinuity 
propagation in second order hydrodynamics and Israel-Stewart theory 



Vdisc = \ hr, — -i — 7 ~ 0.618546 . (42) 
V 2(2 -In 2) v ' 

Next, we solve for the shock solution for speeds that are close to the speed of sound. To 
this end we have to use the solution of first order hydrodynamics (25) and expand to second 
order in u m , the difference between the actual asymptotic speed and the speed of sound. 
Using (23) we now solve 

rptx rptx rpxx rpxX (43) 

Again the temperature field is given by (20) and the second term in the expansion of the 
velocity must satisfy 

Su' 2 = ^Uui (V26 Ul {l - (log 4 - 2)<?«i)+ 

(44) 

+ 65u 2 + V2(\ogA - 7)) + Ay/2 ' 



where, as before, the derivatives are with respect to £ = ^-u^Tix. The solution is given by 



MO = I 



In cosh, c ( £ i 

4>/2(l — In 2) ^ + 5^ tanh 2 £ + tanh£ + — 

cosh £ V c °sh £ / 



(45) 



This solution agrees with the one derived in gravity in section (3), using the prescription of 
[50]. It also agrees with the linear analysis carried out above. Indeed, we can determine q 
through 

iq± = lim (46) 
x^±oo du/dx 

and compare them to the values following from (40). Since this derivation is identical to the 
one we use in section (3.1) we omit it here. 

Finally, we examine the entropy production for this solution. The entropy current has 
been derived in [54, 55, 52]. For second order hydrodynamics the current and the entropy 
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Figure 5: Entropy production as determined by the divergence of the first order entropy 
current s M = Aivqu 11 is plotted versus £ = 4 "°°" T . For the first order solution (25) (higher 
curve) the production is symmetric with respect to the front. For the second order solution 
(45) (lower curve) the entropy production is slightly displaced towards the supersonic side 
of the jet. The normalization of the vertical axis is such that the difference As x for the two 
sides for the first order solution is 1. Notice that the two curves do not have to integrate 
to the same number, since the asymptotic velocity and temperature differ for the first and 
second-order solutions, cf. Eqs. (25) and (45). 

production are given by 

s» = 4vr^ - T -jlo™o KV u» , (47) 

V = ■ ( 48 ) 

It is easy to check that the solution (45) satisfies (47)- (48). Interestingly, as seen in figure 
(5), the entropy production is larger on the supersonic side of the wave. 

Similarly we can examine the asymptotic tail of shock waves in the Israel-Stewart theory 
[56, 57]. This is a theory originally proposed to cure the instantaneous propagation of 
discontinuities in first order relativistic hydrodynamics. A new tensor is introduced that 
parametrizes the departure from the ideal fluid: 

T" v = p {rf v + 4uV) + ir» v . (49) 
The tensor n^" is connected to the velocity and temperature fields by 

7T^ + T^VyJX^ = -2 W ^ + T W (a;^ + C^TT^) , (50) 
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where T>\ is the so called conformal derivative 3 and T n (or alternatively rj) is a parameter 
the value of which has to be determined from microscopies. One often uses the rescaled, 
dimensionless parameters fj and f n defined by 



V = v(nTf, r n = ^. (51) 
For the M = 4 superconformal plasma 

2 -In 2 , ro , 

V = l, *W = — ^— . (52) 

Alternatively, one may initially leave these parameters undetermined and then choose them 
to fit specific quantities. In the case of a shock wave = 0. In order to examine the 
asymptotic falloff in a linearized theory we perturb the asymptotic values of u, T with (37) 
and, in addition, one component of the ir^ u tensor with 



7l°°(x) 



5ir 00 e iqx . (53) 



The resulting three by three system has a solution only if q is given by 



iq 3v 2 - 1 VT 



v 



2 



irT v rj + 7V (1 — 3t> 2 ) 



(54) 



Notice that there is a pole, which for the M = 4 values (52) is located in v ~ .92. On the 
other hand, the fully microscopic calculation based on the gravity dual (and described in 
section (4)) shows that q remains finite at all v < 1. This may lead one to choose fj, f n in such 
a way the pole in (54) is located at v = 1. However, even with this choice the Israel- Stewart 
theory fails to capture the asymptotic behavior of q at v — > 1. Equation (54) predicts that 
q increases linearly with 7 = ^tj=j , whereas the linearized gravity analysis of section (4) 
predicts q oc 7 1 / 2 . 



2.4 Effective Hydrodynamics 

Effective hydrodynamics is the approach where one attempts to model the effect of higher- 
order terms in the gradient expansion of T^ u with terms that are high in derivatives but 
linear in velocity. Such an approach has been taken up in Ref. [58] and may seem an ideal 

3 For a definition of T>\, to' 11 ' and a comparison between the Israel-Stewart theory and second order hydro- 
dynamics for conformal plasmas one can consult [54] 



11 



way to encode results of a linearized theory. Effective hydrodynamics for a conformal theory 
in flat space can be summarized by writing the stress tensor as 



T>™ = prf" + (p + p) u»u v + IF" , (55) 

W v = -2 J dt' J d 3 x'D(x -x',t- t')a^{x', t') , (56) 

D(x,t) = J dcod 3 ke~ iuj+ikx r](u, k 2 ) . (57) 

The effective viscosity t]{oj, k 2 ) is taken to be a function of oo, k 2 such that it correctly 
reproduces the location of the poles of the scalar, shear and sound modes up to the desired 
order. For the M = 4 plasma, it has been calculated up to the fifth order in [58] and found 
to be 

77 = 770(1 + ir) ,iu + m,2U 2 + if] 2) \0Jk 2 + ir] - 3 u 3 + + r} % iu 2 k 2 + r] 0A u 4 H ) , (58) 

P+P 010 1 

% = , ^0,l= 2 - ln2 ; V2,0 = ~^, 

77o, 2 ^-1.379, 772,1 ^-2.275, 770,3 ~ -0.082 , ( 59 ) 

7/4,0 « 0.565 , 770,4 ~ 2.9 , 772,2 ~ 1-1 , 

where we have not given the uncertainties of each coefficient. One can attempt to resum 
this fifth-order expression into a rational expression with one or two poles [58] . 

To compute the asymptotic tails of the shock in effective hydrodynamics, we consider 
again perturbations of the type (37) in the rest frame of the shock. The result is given by 

H 3?j 2 - 1 r 7/0 

7rT = ~^ VY ^ ^T^ -A) ' (60) 

where v is the speed of the fluid relative to the shock. The actual curve q = q(v) is determined 
by linearized gravity and can only be found numerically; the result is shown in Fig. (6). A 
simple ansatz for the effective viscosity does not reproduce this curve very well. For example 
an effective viscosity with one or two poles will always give q — > as v — » 1, in contrast 
to the behavior following from linearized gravity. On the other hand, as seen in Fig. (6), 
the expansion of Ref. [58] gives a reasonable approximation for q on the subsonic side of the 
shock. 

We can take the idea of effective hydrodynamics one step further by simply encoding 
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our numerical curve into an expression for T^", so that agreement with linearized gravity 
is perfect by construction. In general, any hydro dynamic approximation amounts to recon- 
structing T^ v from its timelike eigenvector (cf. Eq. (16)). Consider the fluid at rest with 
a sound wave of small amplitude e propagating along x with momentum k and (in general 
complex) frequency u(k). From conservation of T^ u and the traceless condition we find that 
the Fourier components of 5T^ V are given by 

2 / 2 2 

5T tt = e, ST tx = %e, ST XX = 5T yy = 5T ZZ = ~" e . (61) 

k ' k 2k 2 

It is now a simple matter to compute the timelike eigenvector of T 7 ^, identify 5u = as 
the x component of the four velocity and write 5T^ U as 

ST n„ = 6T t» + 5T »v ^ 

where ST^ is the variation of the ideal fluid energy momentum tensor = (vrTo) 4 {rf v + Au^u u ) 
and is an extra contribution given by 



5T W = ^-{^ 2 -k 2 )T*5u, (63) 



STfi = 5T^ = -^{^-k 2 )T^u. (64) 

In these expressions u should be understood as a function of k obtained by solving numer- 
ically the gravity equations for the sound wave. Notice that the function 3 ^~ fc2 is regular 
for k 0, so Tff. is a well defined function of 8u. It is non-local since it involves an in- 
finite number of derivatives. Nevertheless, in the linear approximation we can work with 
this energy momentum tensor that reproduces exactly the sound pole and the asymptotic 
behavior of the shock wave far from the shock. Later, in the numerical section we give an 
approximate result for the function u(k) that could be used, if so desired, to further simplify 
and approximate ST^. Such effective hydrodynamics is still not sufficient to describe the 
center of the shock, where the linearized approximation is not applicable, but at least it 
summarizes all the information we were able to extract from gravity without attempting to 
find the full numerical solution to the Einstein equations in the bulk. 
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Figure 6: The solid (blue) line represents the numerical results discussed in section (4). The 
(purple) dashed line is the first order hydrodynamics result and the (green) dashed line is the 
second order result. The (blue) dotted line is the Israel- Stewart result and the dotted-dashed 
line is the result for the effective hydrodynamics theory of section (2.4). None of the theories 
captures the v — > 1 asymptotics of the numerical result. 

3 Shock waves in the gravity-hydrodynamics correspon- 
dence 

In a strong shock wave, the region near the shock is beyond the hydrodynamic approximation 
and can only be described as a jump in the hydrodynamic quantities. In the case of the 
strongly coupled M = 4 plasma, hydrodynamics ceases to be valid at distances shorter 
than the inverse temperature. However, at those distances the bulk description in terms 
of gravity does not break down suggesting that gravity can resolve the shock waves and 
provide a smooth description for them. However the velocity and temperature are not well 
defined in the region of the shock, so the best characterization is in terms of the energy 
density, namely T tt (x). In this section we first study such function in the case of weak or 
hydrodynamic shocks. In that case we can reproduce the results of the previous section and 
obtain the dual metric (within the hydrodynamic approximation). Afterwards we consider 
strong shocks and, by using linearized gravity, obtain the width of the shock, as determined 
by the exponential tails on both sides. This is detailed information that can only be obtained 
from a microscopic theory of the system. In principle, one would like to go further and obtain 
T tt (x) by solving the Einstein equations in the bulk numerically, but such a calculation is 
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beyond the scope of this paper. 



3.1 Weak shocks: an explicit solution in the fluid-gravity corre- 
spondence 

The four conservation equations d^T^ = are not enough to determine the nine independent 
components of T M1/ in the boundary theory. Hydrodynamics amounts to a restriction on 
T^ v by providing an expression for it in terms of four variables, the velocity ^=1,2,3 and 
temperature T, which are then determined from the conservation equations. The expression 
for T^ v is given as a derivative expansion, and its precise form can only be determined from 
a microscopic theory of the system. 

From the dual gravitational point of view, the energy momentum tensor sets the boundary 
conditions for an asymptotically AdS metric and the conservation equations d^T^ = are 
necessary consistency conditions for the existence of a solution to the Einstein equations 
with such boundary data. As recently explained by Bhattacharyya et al. [50] those metrics 
generically have naked singularities. The condition of the metric being non-singular imposes 
a restriction to T MV , which is a counterpart of the restriction seen in the hydrodynamic 
construction. In fact, when this analysis is done in a derivative expansion, as shown in [50], 
it provides a microscopic derivation of the hydrodynamic equations for the strongly coupled 
M = 4 plasma. 

In this section we use the BHMR construction [50] to obtain the metric dual to the 
shock waves in the hydrodynamic regime up to terms which are third order in the derivative 
expansion. 

Let us start by summarizing the procedure as adapted to our particular problem. The 
starting point is the boosted black brane in Eddington-Filkenstein coordinates: 

(ttT(x)) 4 

dsf 0) = -lu^dx^dr + \^u lL u v dx l "dx v + r 2 rj^dx^dx u , (65) 

where 

= (n°(x),n(x),0,0), u°(x) = y/l + u(x) 2 . (66) 
Since u(x) and T(x) are not constant this metric does not solve the Einstein equations 

G M n = Rmn + ^9mn = , (67) 
where M,N — r, t, x, y, z denote 5-dimensional indices (whereas v = t, x, y, z denote four- 
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dimensional indices). To find a solution we introduce a formal parameter e and expand 

T(x) = T ( 0)(x) + er (1) (x) + e 2 T (2) (x)+e 3 T (3) (x) + O(e 4 ) . (69) 

At the same time, the metric is corrected by adding an expansion 

ds 2 = ds 2 {0) + eds 2 (l) + e 2 ds 2 2) + £>(e 3 ) . (70) 

The tensor Gmn is expanded in powers of e with the rule that each ^-derivative counts as one 
extra power of e. The equations are solved order by order. For that purpose it is convenient 
to introduce the vector 

u» = (u{x),u°{x),0,0) , (71) 
orthogonal to u^. In this way the metric corrections can be parametrized as 

^ s (n) = ( s i + s^)u fl u u dx^dx u + s^T] tJ-u dx fJ 'dx 1 ' + 2s^u tl dx^dr + Sy* > u tl u l/ dx fJ 'dx 1 ' 

+4 n) [u y u v dx i "dx v - l -{dy 2 + dz 2 )^j , (72) 

where the functions s[ n 2 3 ,sy,sr describe scalar, vector and tensor perturbations classified 
according to the local SO(3) group that leaves w M invariant. Notice that each ds 2 n ^ should 
in turn be expanded using (68). Following [50] we make the gauge choice gffl = 0, g$ ~ 
and g^guj} = for all n > 0. In that case it is convenient to parametrize the fluctuations as 

s$ n) = ^JfcWfor), 4 n) = r 2 /iM(x,r), = §/^(x,r), 

Sv = ^j^{x,r), 4 n) = rV n )(x,r) , (73) 

In order to solve the equations order by order it is convenient to decompose Gmn into its 
scalar, vector and tensor parts which decouples the equations. The procedure is in principle 
straightforward, and we proceed to describe the results. 

Order 1 . The first equations that we obtain are u' = 0, T/ x = implying that 

M (0 )(x) = M( ), T (0) (x) = T(o) , (74) 
namely they are constant functions. In that case the zero order metric is an exact 
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solution and there is no first order correction to the metric: 



/l (i) =A ,(i) =J (i) =a (D =0 . (75) 
The temperature however is corrected to 

= -^T (0 ,«(i) • (76) 

Order 2 . The first equation we find is 

M ' (1) (2^ 0) -l)=0. (77) 

We can take itm constant which leads us to a trivial solution, or otherwise we require 
U(o) = implying that the zero order solution describes a fluid moving at the speed of 
sound, which is the appropriate starting point to describe shocks in the hydrodynamic 
approximation. The other equations give: 

M 2 ) = 0, fcW * V j^ = -^r\' a W = ^-F 1 (-? r -) , (78) 

3 (1) ' v/3 37rT {0) \nT (0) J 1 ; 

where 

F 1 (y) = ln^ (1 + y ^ 1 + y) ^ -2arctany + 7r. (79) 

Notice that at this order um is undetermined. This is a particular property of our 
solution that requires us to go to higher orders to obtain the metric. 

Order 3 . The first equation we obtain is 

wji) = 3 7rT (o)«(i)«(i) , (80) 

which allows us to solve for urn as we do further below. The components of the metric 
are corrected by 



J {3) = — ^%) + ^W(i)(^ (0 )) 3 F 2 



r 



(0) 



,(3) = (_ 

37rT (0 ) V' 
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with 



F2(y) = ^(y 4 -l)(2arct a ny + ln(^±^^-iv + y 3 + 2/ 2 -§, (82) 



F'M = ^^{2(l-y 3 )(arctan 2 /-^-ln(l + y)) + (l + y 3 )ln(l + y 2 ) 



-% 3 lny + 7-21n2-^ : ±|l- T ^--V) . (83) 

i + y 2 l + y J 

We give the derivative of _F 3 since it is the function that enters in subsequent calcula- 
tions. It can be integrated explicitly in terms of dilogarithms but the expression is not 
very illuminating. The temperature is given by 

nT {2 ) = -^7rT ( o)U (2 ) + ■ (84) 

We want to write the metric up to order e 2 which requires computing it( 2 ). Since it is 
undetermined at this order we continue the expansion. 

Order 4 . At this order we only look for the equation determining uq). However we need 
to include and keep track of the terms W( 3 ), T( 3 ) , etc. to be sure that they do not 
appear in such equation. What we get is 

u" {2) - -ttT (0) (u (1) u (2) )' = ^-(7- 21n2)u' (1) 2 - -^^0)^^(1 + 4^2) . (85) 

The last equation, together with (80), can be easily solved to obtain, at this order 

u = W( ) + U(i) + U(2) , (86) 

U( ) = (87) 



U(i) = -Mootanh^, 

In cosh. £ ( £ 

±y/2(l - In 2) + 5^ tanh 2 £ + tanh f + 



u 2 

"oo 



«P0 = 6 



cosh" 1 £ V ' cosh 2 £ 
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where 

e = — f-^x, (90) 

and we set the formal parameter e = 1. It is interesting to note that the equation for Mm 
determines that the fluid reaches the shock wave supersonically and leaves subsonically. In 
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other words, we are not free to exchange the subsonic and supersonic sides of the shock. The 
reason is that are choosing gravity solutions which are regular in the infalling Eddington- 
Filkenstein coordinates as appropriate for a black hole. The constant is arbitrary and 
determines the amplitude of the shock, namely, the asymptotic value of the velocity. For 
consistency of the approximation we require ii m <C -j=. In fact, plays the role of the 
small parameter, as can be seen from the fact that the velocity depends on x through £ and 
so each x derivative brings in an extra power of Uoo. The behavior at infinity is given by 

u(x-+±oo) ~ ^Tu 00 ±2u 00 e**\l + ~u 00 (4V2-4]n2±5V2)^ (91) 

We write the correction in the exponential form for an easier comparison with the hydrody- 
namic result. For the temperature we have 

T(x -> ±00) = T( ) ± — T (0 )Moo • (93) 

For the velocity we have 

1 1 2^2 14 2 
v(x ->■ 00) ~ —= + ov + = — -Uoo H F^oo , (94) 

1 1 2^2 2 2 

v{x -> -00) ~ —= + ot>_ = —= H ^Moo -j=u . (95) 

\/3 \/3 3^3 3^ 

The first check is that the condition (6) 

v{-oo)v{+oo) = ~ + 0(i4) , (96) 

is satisfied to the considered order. From eq. (92), the exponential tail is given by 

v(x -> ±00) = -= + £i;± + c ± e i9±x , (97) 
v3 

for some constants c± and the penetration depth determined by 



iq± 



2y/65v± + 6y/2Svl(l - In 2) , (98) 



7TT ± 

in complete agreement with the hydrodynamic calculation. This is not surprising since the 
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hydro dynamic equations arise from gravity. The calculation in this section, however, allows 
us to compute in addition the dual metric, which to this order is given by 

ds 2 = -lu^dx^dr + -^((vrT) 4 + k)u ll u v dx'*da? + r 2 r]^ u dx^dx v + \ju p u v dx il dx v (99) 
+ar 2 (u^ 1 u v dx ll dx u - ^{dy 2 + dz 2 )) , (100) 



where = (\/T+ u 2 , u, 0, 0), = (u, y/l + u 2 , 0, 0) and u should be expanded as u 
M( ) + + U(2) using the U( n ) computed in (89). The temperature is also expanded as 



T = T(0) ~ 1T T(0)U(1) ~ ^Y T ^ u w + > ( 101 ) 



whereas the other functions entering the metric are given by 

k = ^r 3 (n' (1) + n' (2) )-^r 2 ^), (102) 

2 4\/2 / r \ 

i = -^^ 3 Ki) + ^2)) + ^(i)Ki)(^ (0 )) 3 F 2 ^-^-J, (103) 



u '(i)+ u '(2) ( r \ Ay/2 
37rT ( o) 1 \ttT {0) J 97rT {0) "^"^ J V.^(o) 



« = ^fe^* bb " db"««<l>* bb • ( 104 ) 



with the -^1,2,3 as defined above. The expansion parameter is the strength of the shock as 
determined by the constant appearing in u< n y 



4 Strong shocks: the linearized gravity approximation 

Strong shocks are characterized by large gradients of velocity and temperature and cannot 
be studied within hydrodynamics: we can say that hydrodynamics does not resolve their 
profiles. However, the gauge-gravity correspondence is not limited to small gradients, and 
so the gravity side of it should contain information about strong shocks as well. We now 
discuss how that information can be extracted. 

In principle, we expect that there are exact solutions to the 5-dimensional Einstein equa- 
tions, and those solutions describe strong shocks exactly. This belief is based on the obser- 
vation that the asymptotic values of and T on the far left and far right of the shock have 
vanishing gradients and are therefore well reproduced even by the ideal hydro (see Sec. 2). 
On the gravity side, to each of these asymptotics, there corresponds a 5-dimensional AdS 
black brane, suitably boosted and with a suitable value of the temperature. Then, there 
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must be a 5d solution describing a stationary wave that smoothly interpolates between these 
two regions — the gravity dual of a strong shock. 

The exact solution (assuming it exists) described in the preceding paragraph would tell 
us all there is to know about a strong shock, in particular, the profile of the average energy 
density, e(x). So far, however, we have not been able to find any such solution explicitly. In 
this section, we provide partial information about the profile of a strong shock, obtained by 
looking at linearized gravity on the backgrounds corresponding to each of the two asymptotic 
regions (x — > ±00). 

4.1 Equations of linearized gravity 

In linearized gravity, one writes the metric in the form = gj$ + h^ v , where g$ is the 
metric of the AdS black brane, and h^ u is a perturbation, and works to the first order in 
the perturbation. Solutions to the linearized Einstein equations are known as quasinormal 
modes. In this paper, we consider solutions that depend only on the coordinate (x) along 
the direction in which the shock propagates and, possibly, time. In this section we adopt 
the convention ttT = 1, the temperature dependence can be recovered by multiplying the 
modes by irT. Due to the translational invariance along the brane directions, we can search 
for these solutions in the form 

h^(t,x,r)=r 2 H^(r)e- M ^. 

We adopt the convention that u and q refer to the boosted frame, moving at the speed of the 
shock, and their primed counterparts, u' and q', to the unboosted frame, connected with the 
fluid. Note that there are actually two such unboosted frames (the fluid on the two sides of 
the shock moves with different velocities), but in the linearized approximation the two sides 
are disconnected and can be considered separately. 

Linearized gravity has by now become a familiar tool in studies of the kinetics of the 
strongly coupled M = 4 plasma but in a setting that is typically different from ours. In many 
cases (as, for example, in the computation of the viscosity [51]) one considers relaxation of an 
initial perturbation. Then, one picks a real wavenumber q and looks for the corresponding 
(complex) frequencies. Here, in contrast, we are interested in propagation of a boundary 
disturbance, that is, in how far a perturbation with a given frequency extends into the 
plasma on either side of the shock. For this, we pick a real u and look for the corresponding 
(complex) q. Specifically, we will be interested in perturbations with u = 0, as we expect 
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these to describe the behavior of the average energy density of a shock wave sufficiently far 
away from it. 4 

A classification of the quasinormal modes of an (unboosted) AdS black brane has been 
given in [59]; a comprehensive recent review of quasinormal modes is [60] . According to 
that classification, metric fluctuations group into several channels, corresponding to different 
gauge-invariant combinations of the components of H^ u . Here, we are interested in the sound 
channel. The corresponding quasinormal modes satisfy the equation 



Z" + P(u)Z' + Q(u)Z = 0, (105) 

where primes denote derivatives with respect to u = r^/r 2 , 
3uJ ,2 (l + u 2 ) + q' 2 (2u 2 -3m 4 -3) 



P(u) 
Q(u) 



uf(u)[3u' 2 + q' 2 {u 2 -3)] 

4g'V 3u/ 4 + g' 4 (u 4 - Au 2 + 3) + u l2 q' 2 (Au 2 - 6) 



f{u)[3u' 2 + q' 2 (u 2 -3)] 4uf 2 (u)[3uj' 2 + q' 2 (u 2 -3)} 



and f(u) = 1 — u 2 . Remember that in these expressions, u 1 and q' are in units 5 of irT and 
refer to the unboosted frame. They are related to the frequency and wavenumber in the 
boosted frame by the Lorentz transformation 

lo' = uj cosh (3 — q sinh (3 , (106) 
q = —u sinh j3 + q cosh (3 , (107) 

where tanh /3 = v , the speed of the shock. For static perturbations, we set uj = and 
substitute the resulting expressions for u' and q' into Eq. (105), to obtain: 

p( . 3 + 3u 2 - 5 7 V + 3 7 2 ^ 4 nn8 , 

P ( u ) = 77 w 2 2 , (108) 

uj(u){Yu z — 3) 

4t 2, u 2 o /2 m 2 1 

<2( M ) = - 77 w p 9 ^ + g 2 fo/ ■> ) (109) 

/(m)(7V-3) Auf 2 (u) 

where we have used the shorthand 7 = cosh /?. 

The same expressions can be obtained by starting directly in the boosted frame. In this 



4 To be sure, it is not obvious a priori that the shock wave profile is static: there could be instabilities 
in the nonlinear central region that cause oscillating behavior. In the linearized theory, a possible signal of 
such an instability would be the absence of a physically acceptable static solution (due. for instance, to a 
singularity in the equation). We have not found any such signals in our calculations. 

5 And so are twice as large as their counterparts in [59]; hence an extra overall \ in the second term in Q. 
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case, the unperturbed metric is 



ds = r 2 rj liv dx^dx y H — \{dt cosh (3 — dx sinh (3) 2 1 



r 2 r 2 (l — r 4 /r 4 ) ' 

and the perturbation reads 

ds\ = r 2 [H 00 dt 2 + H n dx 2 + 2H 01 dtdx + #(A/ 2 + dz 2 )] e iqx 

(all if^,, are functions of r only). The relevant gauge-invariant combination (at u = 0) is 



r, 



Z(r) = #oo(r) + ( 1 + '\l 2 ) # (r) , 



r 1 



and satisfies Eq. (105) with the coefficient functions given by Eqs. (108) and (109). 

For computation of properties of the plasma via the gauge-gravity correspondence, we 
only need to consider Eq. (105) in the region outside the horizon, r < r < oo, which we will 
refer to as the physical region. In terms of the variable u, it corresponds to < u < 1. A 
noteworthy property of the coefficient functions (108) and (109) is that, at sufficiently large 
boost velocities, 

cosh/3 > cosh/3 cr = \/3 , (HI) 



P and Q both have poles inside the physical region, at 

cosh/3 



u = U\ 



(112) 



From the outset, we might have anticipated that we would need to impose boundary con- 
ditions at the boundaries of the physical region, u — > and u — > 1, but not at any interior 
point. We therefore need to explore the nature of the singularity at u = u\ in more detail. 

Let us search for solutions near u = u\ in the form Z ~ (u — u\) s . For the exponent s, 
we find two roots, 

s = or 3. (113) 

Fuchs's theorem [61] guarantees that the larger root correspond to a regular solution, expand- 
able in powers of w = u — U\ as follows: Z = w 3 + 0(w 4 ). As for the solution corresponding 
to the smaller root, in general, we expect it to have the form 

Z{u) = c + ciw + c 2 w 2 + c 3 w 3 + c' 3 w 3 \nw + . . . (H4) 
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The recursion equation for the coefficients c n , which is obtained by substituting Eq. (114) 
in Eq. (105), degenerates at the order (and only at the order) at which the second solution 
appears, in our case the order w 3 . Whether or nor the logarithmic term in (114) is nonzero 
then depends on the precise values of the coefficients of all the terms up to 0(w) order in 
the expansions of P(u) and Q(u). As it turns out, there is a curious cancellation among 
these terms, such that c' 3 = for all values of 7 and q. We conclude that both solutions are 
regular at u — u±, and a boundary condition there is not required. 

4.2 Boundary conditions. Irreversibility 

The choice of boundary conditions for quasinormal modes that is suitable for applications 
of the gauge-gravity duality to kinetic theory has been discussed in the literature (see, for 
example, Ref. [51]), and we do not deviate from it here. Our computation, however, re- 
quires an analytical continuation of these boundary conditions, which is the subject of this 
subsection. 

At u — > (the boundary of the AdS space) we use the standard 

Z(0) = 0. (115) 

At u — > 1 (the near-horizon region), we first consider real u' and pick, as usual, the wave 
infalling with respect to the black brane: 

Z(m->1)~(1-m)"^'/ 4 . (116) 

For the present problem, having to do with propagation of a perturbation in space, rather 
than in time, we need to analytically continue this expression to complex u' given by the 
Lorentz transformation (106) (u/ is complex because so is q). In particular, for the static 
case (oj — 0), we have u' = — gsinh/3 and thus 

Z(u ->■ 1) ~ (1 - u )¥i^P . (117) 

This choice is equivalent to choosing the solution that is regular in infalling Eddington- 
Filkenstein coordinates as done in section 3.1. There is an exceptional case (3 = (3 cr , the 
critical value given by Eq. (111). For this value of 0, the analytical continuation to w = 
causes confluence of the singularities at it = 1 and u = u\, which modifies the asymptotic 
behavior near u — 1. We consider this case separately later in this subsection. 
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Recall that Z(u) corresponds, via the gauge-gravity duality, to a perturbation in the 
plasma that depends on x as exp(iqx). On physical grounds, we expect that perturbations 
corresponding to the tails of a shock at x — > ±00 decay away from the shock. This means 
that we must pick q with a positive (negative) Im q for the fluid at positive (negative) x. 
Recall also that positive (negative) x correspond to the subsonic (supersonic) side of the 
shock. Thus, according to Eq. (117), Ziu) is regular at the horizon on the supersonic side 
but singular on the subsonic one. 

Formulating a boundary problem for the regular case presents no difficulty: the second 
solution to Eq. (105) diverges at u — > 1, and the boundary condition (117) selects the one 
that does not. The singular case (Img > 0) is a bit trickier: we wish to retain the divergent 
solution and reject the convergent one. To achieve that, we peel off the singular part, as 
follows: 

Z{u) = (1 -u) s ^(u), (118) 

where 

1 1 

s = ~^ LUJ ' = -igsinh/3 (119) 
(Res < 0), and demand that ty(u) is analytic at u = 1. This works whenever 

2s 7^ integer . (120) 

Indeed, the regular solution behaves as (1 — u)~ s , and the corresponding ^ as (1 — u)~ 2s . 
Provided the inequality (120) is satisfied, this is not analytic and is rejected by our boundary 
condition. 

Note that the inequality (120) is sufficient but not necessary for the singular boundary 
problem to make sense. Suppose (120) is not satisfied for some values of q and /3, but both 
solutions for ^ are regular at u — 1. We consider such a q to be an eigenvalue of our problem 
(at that particular /?), because we can always form a linear combination of the two regular 
solutions that satisfies the second boundary condition (115). On the other hand, it is a priori 
possible that (120) fails in such a way that one of the solutions for Z contains a logarithm 
of 1 — u. In this case, we truly have no recourse; indeed, the singular part cannot even be 
peeled off as in (118). Interestingly, the condition (120) never breaks down for the "main" 
branch of q(v), as defined below. 

We anticipate that there is more than one eigenvalue of q for each value of the shock's 
speed v. We refer to these as different branches and denote them as q n (v). Let us mention 
some of the properties of these eigenvalues for the case 00 — 0. 
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First, setting q = —in makes all the coefficients in the equation (105) real and turns the 
condition (117) real as well. We conclude that the eigenvalues q n (v) are all purely imaginary, 
and s all purely real. (This is not the case at u 7^ 0.) 

Second, while the functions (108) and (109) do not depend on the sign of q, the condition 
(117) does. Hence, the set of the eigenvalues q n (v) at a nonzero v is not symmetric about 
q = 0. The reflection q — > —q, without changing the direction of the shock's velocity, is 
equivalent to reflection of both space and time: x — > —x and t — > —t. In particular, it 
exchanges the subsonic and supersonic sides of the shock. The absence of symmetry under 
q — > —q corresponds to the condition (already noted in Sec. (2.2)) that the fluid must be 
supersonic in front of the shock and subsonic behind it, and never vice versa — a condition 
that reflects, ultimately, the second law of thermodynamics. In calculations on the gravity 
side, the source of this irreversibility is the choice of the infalling wave in Eq. (116). 

The main branch. The branch qo{v), for which exp(iqx) decays away from the shock the 
slowest, will be referred to as the "main" branch and often denoted simply as q(v). This 
is the branch that crosses zero at v = 1/ y/3 and is the only one seen in the small-gradient 
(hydrodynamic) approximation discussed in Sec. 2. 

The exceptional case. The preceding discussion of the boundary conditions does not 
apply to the exceptional case u = 0, /3 = /3 cr , when the singularities at u — 1 and u = u\ 
coincide. This case needs to be considered separately. At uj = and (3 = f3 CT , Eq. (105) 
becomes 




(121) 



Solutions near u = 1 are of the form Z(u) ~ (1 — u) s with 



s = 



l±i 



2V2' 



(122) 



Consider the substitution 



Z{u) = (1 - u 2 ) 1/2 $(n) . 



(123) 



Eq. (121) becomes 




Note that for q 2 = —2 the coefficients in Eq. (121) are all regular at u = 1. Hence, q = ±i\/2 
are eigenvalues of the boundary problem. Since /3 = /3 cr is supersonic, only 



q = -iV2 



(125) 
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is physical; it lies on the main branch. Curiously, although the behavior at u = 1 prescribed 
by Eq. (122) is in general different from that prescribed by Eq. (117), for q = —i\/2 (and 
s = 1/2) they coincide. As a result, the curve q(v) corresponding to the main branch is 
continuous at v = v CT . 

For branches above the main branch, iq n (v CT ) is large, so that only one of the solutions 
to Eq. (121) is regular at u — 1. The boundary condition is to choose the regular solution, 
which corresponds to choosing the plus sign in Eq. (122). This is equivalent to using 

Z(u) = (1 - u) i<?/V %(«) (126) 

[cf. Eq. (117)] but allowing ty(u) to vanish (linearly) at u — 1. Indeed, as v goes through 
v cr , \I/(1) goes continuously through zero. Thus, these other branches are also continuous at 



4.3 Quasinormal modes for special values of v 

4.3.1 v = (fluid at rest) 

The case of a plasma at rest has already been studied. In particular, as argued for example 
in [62], the eigenvalue q is given by the lowest glueball mass in Witten's QCD3 construction 
[63]. The reason is that, in the Euclidean space, both finite temperature M = 4 and QCD 3 
are dual to the same AdS black hole. For the channel we are considering, the glueball mass 
was computed in [64] giving iq = 2.3361 in perfect agreement with our numerical results. 
This provides a nice check of the calculation although we should point out that, in the frame 
where the shock wave is at rest, the velocity of the fluid is always v > | so v — is not 
directly relevant to our problem. 

4.3.2 v = l/y/3 (the speed of sound) 

This is the limit when the strength of the shock (as measured by changes in various quantities 
between the left and right of the shock) vanishes. As we have seen in section (2), in this 
limit q = 0. The corresponding solution to Eq. (105) is 

Z{u) = u 2 . 
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4.3.3 v = a/2/3 (the singular point) 

In this case the pole corresponding to the horizon u = 1 merges with the one at u — u\ 
since v = a/2/3 implies u — u\. Interestingly in this case we can find the exact eigenvalue 
iq = a/2 and the eigenfunction is given in terms of a hypergeometric function: 

= mv/T^ l — J F 2>1 ( — , — ^— ; 1 + i; 1 + -J . (127) 
With this definition we find 

Z{u) = -Re{y{u)) + Cilm(y(u)) , (128) 

where C\ is a constant that we evaluate numerically to be C\ = 0.38898 from the boundary 
conditions. 



4.3.4 v — > 1 (ultrarelativistic limit) 

Numerical solution (described in the next subsection) shows that the values of s, Eq. (119), 
for the physical branches become large in the limit v — > 1, and the maxima of the eigenfunc- 
tions scale towards u — 0. This suggests that we can obtain an equation applicable in the 
ultrarelativistic limit by neglecting u in comparison with unity in the coefficient functions 
(108) and (109). We obtain 

™ - • (129) 

Q(u) = - ^_ + f £?zl. (130) 

/ -f z u z — 3 4u 

These expressions suggest further that we define a new variable x, as follows: 

u 2 1 

x = = -u 2 cosh 2 (3 , (131) 
uf 3 

and take the formal limit cosh /3 — > oo while keeping x fixed. Eq. (105) becomes 
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where primes now denote derivatives with respect to x, and 



2 2 

p = q m — 




(133) 



The change of variables (131) and the limit cosh/3 — > oo map the physical region < u < 1 
to < x < oo, with the Dirichlet boundary conditions at both ends. As before, the point 
x = 1 (formerly u = u±) is a singular point of the equation but not of either of the two 
linearly independent solutions. Thus, in numerical integrations we can circumvent this point 
by first displacing it into the complex plane, i.e., replacing 1 — x in with 1 — x + ie in Eq. (132) 
(the sign of e does not matter), and then taking the limit of the solution at e — > 0. 

It follows from Eq. (132) that all solutions must have extrema (maxima, if we agree to 
choose the overall sign of Z in a certain way) at x — 1. The solutions vanish linearly at 
x = and, provided Im p < 0, exponentially at x — > oo: 



In Fig. 7, we plot the eigenfunctions corresponding to the smallest two values of |Im p\. 
According to Eq. (133), these determine the asymptotics of the main (lowest) and the next 
lowest branches of q n (v) in the ultrarelativistic limit v — > 1. Numerically, we find 



Note that numerical solution is needed only to determine the coefficients in these formulas: 
the scaling with 7 follows directly from Eq. (133) and the fact that Eq. (132) contains no 
parameters. 

4.4 Numerical results 

Apart from the very few values of v (discussed earlier) for which we have found analytical 
solutions to Eq. (105), we have resorted to solving this equation numerically. We have used 
two numerical methods: (i) the shooting method and (ii) the series expansion. Where their 
domains of applicability overlap, these methods have produced equivalent results. 

In the shooting method, we peel off the non-analytic part as in (118) and set up an initial 




(134) 



iqi{v 



iq {v 



1.895^7, 
5.424^. 
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Figure 7: The wavefunctions (arbitrarily normalized) corresponding to the ground state 
(solid line) and the first excited state (dashed line) of Eq. (132). 

value problem for ty(u) at u close to 1. From (105), the expansion of ^f(u) near u — 1 is 

Z(u -> 1) = (1 - w) s [1 + A(l - u) + . . .] , (135) 
where s is given by (119) and 

8(2s + 1) 
The initial value problem is 

= 1, (137) 
$'(1-5) = -A. (138) 

We can then adjust q (on which both q' and u' depend) so that the boundary condition 
(115) at the other end is satisfied. The limit 8 — > is expected to be smooth whenever the 
boundary conditions (137)-(138) are sufficient to reject the second solution. This is always 
the case for the regular problem (Res > 0), but not for the singular one (Res < 0). In the 



/2 



As + 



16^(25 
2q> 2 - 3cu 



(136) 
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latter case, (137)-(138) are sufficient only if 

Res>-±, (139) 

which is a stronger condition than (120). 

Another caveat is that we need to develop a way for circumventing the singular point 
u = Mi, in the case when the shock velocity (relative to the fluid) exceeds the critical 
value given by (111), and the singularity moves into the physical region < u < 1. Even 
though, as we have seen, solutions to Eq. (105) are always regular at u — Ui, the singularity 
in the coefficient functions precludes passing through this point by means of a numerical 
integration. The approach we adopt here is to consider solutions that are not exactly static 
in the boosted frame but oscillate with a small (real) frequency oj. A nonzero u displaces the 
singularity into the complex plane, so that the equation can be integrated numerically. The 
eigenvalues and eigenf unctions at u = can then be obtained as limits of those at u ^ as 
uj — > 0. The absence of singularity in the solutions guarantees that these limits are smooth. 

In the series expansion method, one develops two power series expansions, one near u — 0, 
the other near u — 1, starting with the terms prescribed by the boundary conditions (115) 
and (117) (after peeling off the non-analytic behavior at u — > 1 as in (118)). The logarithmic 
derivatives of these two expansions are then matched at an intermediate point of the interval 
< u < 1. One expects that, if the expansions near the endpoints are taken to sufficiently 
high orders, the results will be insensitive to the precise value of u at which the matching 
occurs. This method does not require any special device to circumvent the singular point 
u — u\. Indeed, that can be verified by developing a third series around the point u = u\ 
and then matching the logarithmic derivatives with the two series developed around u = 0, 
u — 1. The results are indistinguishable numerically. 

In Fig. (8) we show several branches of q n (v ) obtained by these methods. It is interesting 
to note that a good approximation to the main branch is given by 6 

£-*(§)^(-£)- (140) 

A better approximation can be found by including more parameters in the fit. Including one 
more parameter, the curve 

iq ^1 - Vy/3 (m) 



ttT v 1 'Fv + G 



'Here we restore the dependence on T . 
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-Im(q)/7iT 




at the speed of sound. The dotted sections of the curves denote unphysical values of q — 
those that are discarded as they correspond to perturbations growing exponentially at the 
respective asymptotic infinities. The solid sections denote physical values (corresponding to 
— Im(g) < for subsonic 7 and — Im(g) > for supersonic). The asymptotics for 7 — )■ 00 
for the main branch and the one above it are discussed in section (4.3.4). 



where 

^I(?) 3 V-D>, G= - 3-"(l+ff-*») , (U2) 

gives a better approximation close to the speed of sound and the large 7 asymptotics. 

These approximations can also be used to obtain an approximate function u)'(q') (the 
dispersion law) for the sound waves. Transforming to the unboosted frame and using v = 
— we find that the approximation (140) provides us with the following implicit equation 
for uj'(q'): 

gV 2 -^ 2 ) 3/2 = 16 0) 5 (^ + -^) • (143) 

The other fit can also be used in this way, but the resulting equation is more complicated 
and we omit it here. 



5 Discussion 

In this work we have used the AdS / CFT correspondence to study shock waves propagating 
in a strongly coupled plasma. Shock waves appear quite generically when the motion of 
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the fluid is supersonic and produce dissipation and drag even for zero viscosity. In the 
case of ideal fluids they are associated with surfaces where the velocity and pressure are 
discontinuous. The discontinuity indicates a failure of the hydrodynamic approximation 
and should generically be resolved by a microscopic description of the system which, in 
this context, is provided by the dual gravity description. An exception is the case of weak 
shocks, which propagate close to the speed of sound, where the inclusion of dissipation, 
namely, viscosity, resolves the shock. In this case, the dual metric can be found using an 
expansion in the strength of the shock. On the other hand, strong shocks are beyond the 
hydrodynamic approximation. They can only be resolved by finding the dual gravity solution 
which should be a smooth wave propagating without deformation on the horizon of the black 
hole. Far from the shock, the solution differs slightly from a boosted black hole which allows 
for a perturbative study of the solution. In particular, we have computed, in the rest frame of 
the shock, the exponential tail of the solution, namely, the width or penetration depth of the 
shock. It is a function of the velocity, which we have determined numerically. In particular, 
when the speed of the incoming fluid approaches the speed of light, the penetration depth 
ahead of the shock goes to zero as the inverse square root of the gamma factor, i ~ 7" 2. 
Since the length scale goes to zero, this scaling exponent is an ultraviolet property of the 
theory, as can also be seen from the bulk calculation, where the exponent is determined by 
the properties of the metric near the boundary. It would be of interest for future work to 
establish the value of this exponent for other backgrounds in the context of AdS/CFT or 
perhaps even directly from perturbative gauge theory calculations, for example, in QCD. 
More generically, since shock waves probe microscopic properties of the system, they are 
an ideal tool to study the transition from the microscopic to an effective hydrodynamic 
description. For example, we have shown that, for strong shocks, the dependence of the 
penetration depth on the velocity of the incoming fluid is not correctly reproduced either 
by second order hydrodynamics or by the Israel-Stewart theory. It is possible to encode 
our results into effective linearized hydrodynamics of the type proposed in [58], but with 
effective viscosity given by a numerically determined function of u' and q' . The linearized 
description, however, is valid only far from the shock. It would be interesting to see if an 
improved effective description exists that can correctly capture the main properties of shock 
waves in the nonlinear region. 

Although we have understood several basic properties of shock waves in the context 
of AdS/CFT, there are many interesting questions that we have not addressed here and 
would be interesting to pursue. One important question is if the full solutions dual to shock 
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waves can be found analytically or by numerical methods. They are interesting objects in 
gravity since they correspond to black branes with different asymptotic temperatures on the 
two sides of the wave. Such waves propagate without deformation and generate entropy 
by expanding the area of the horizon. Perhaps they are quite a generic phenomenon not 
restricted only to examples appearing in the context of AdS/CFT. Other, perhaps simpler 
problems to consider are related to the introduction of dynamical quarks by means of probe 
branes [65] in the background of the shock. A shock should appear on the brane giving 
rise to a force on quarks and meson emission from the shock. Finally, the introduction of 
dynamical quarks can also provide a closer point of contact with the quark-gluon plasma 
experiments at RHIC where generation of a Mach cone by a heavy quark propagating in the 
plasma has been recently suggested [66]. For that reason, it would be of great interest to 
understand the conditions under which a moving quark generates a strong shock such as the 
one studied in the present paper. 
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A The equations of motion for the perturbations 



In this appendix we derive the equations of motion for various perturbations of the metric. 
The background metric is given by 

dsl = r^d^dx" + - dx ) 2 + TJ ^r , (144) 

r 2 a/1 - v 2 a/1 - v 2 r 2 (l - r 4 /r 4 ) 

and the perturbation by 

dsj = r 2 e iqx (H 00 dt 2 + H u dx 2 + 2H 01 dtdx + H(dy 2 + dz 2 )) . (145) 
There are seven independent equations coming from 

Gmn = Rmn + 4g MN = . (146) 

We form the linear combination 

£ = A MN G MN , (147) 

where the only non zero entries for the A matrix are the linearly independent equations, 
which are the rr,rt,rx,tt,xx,yy,xt components of Gmn ■ After choosing four entries for 
Amn, namely A rr , A rx , A xx , A yy , one can eliminate Hqi,Hu and their derivatives from C 
Only four constants are needed since two of the equations are first order. After this operation 
£ is a function of only H 00 , H and their derivatives. One cannot use the three remaining 
constants to eliminate one of the functions and it's derivatives. The reason is that only two 
constants are free, the third one can be thought of as an overall rescaling of L and there are 
three coefficients to eliminate, the three factors multiplying H, H', H". A redefinition 



#oo(r) = Z(r) - g{r)H{r) , (148) 

and choosing g(r) such that the coefficient of H" vanishes allows us to write a decoupled 
equation for Z(r). With the choice of 

9(r) = l + r 4^— 2 , (149) 
r 4 1 — v z 
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the final equation for Zir) is 



*"(r) + *(r)i(^=4 + 



r ^ r 4_ r 4 r 4_ 3? ,4( 1 _^ / 

z(r0 f ^ r 4 -r 4 (W) \ 1 J 

^ ^ r 2( r 4 _ r 4)( r 4 _ 3r 4^ _ ^2)) « ( r 4 _ _ „2) ^ 

This equation coincides with the equation for the sound pole [59] when one boosts to the 
frame where the black hole is moving. Now we can trace back the equations and find the 
equations of motion for the rest of the perturbation components. Tracing back the procedure 
to derive the equation for Zir) we find that 



r ( r 4 _ r 4) W r 4 _ 3r 4^ _ ^2) V y ^4 _ r 4)( r 4 _ 3^ _ ^2)) ' 

(151) 

2r 4 (r 4 (l + t; 2 ) - 3r 4 ) , _ , /2 , 4r 4 



# 00 (r) - ^oo(r) r ^ 2(r4 _ ^ _ ^ _ ^ + (r) + 3 (r 4 _ 3r 4 (1 _ ^ 

4r 4 (-2rj + r 4 r 4 (4 - 3t; 2 ) + + t> 4 - 2)) ( } 

+ ^ r (r 4 - r 4 )(r 4 + r 4 (l - v 2 ))(r 4 h - 3r 4 (l - v 2 )) ' 

We can treat the Z(r), Z'{r) terms as a source since Z{r) satisfies a decoupled equation. For 
the last two perturbations H m , Hu it is easier to define a linear combination of them 

H(r) =vH 11 (r) + 2H 01 (r) . (153) 

Treating terms containing H(r), Z(r) as sources H satisfies 

H' {r) -H(r) K + Z\r) ^ " ^ ^ " 2) + ^ " ^ 



(r 4 -r 4 )(r 4 -3r 4 (l -v 2 )) 2 



2ot 4 (3r^ 2 - 2) - 9r 8 (l - t> 2 ) 2 + r 4 r 4 (17 - 20^ 2 + 3^ 4 )) 
(r) r(r 4 - r 4 )(r 4 - 3r 4 (l - t> 2 )) 2 + 

«(j 2 r 6 (l — w 2 ) 

+ Z{r x 



(154) 



-£T(r; 



(r 4_ r 4 )(r 4_ 3r 4 (1 _, ; 2 )) 

v (q 2 r%v 2 - 1) + 2r 4 (r 4 (5 - 3v 2 ) + r 4 (w 2 - 2))) 
r 5 (r 4 — r^)(l — v 2 ) 
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Having determined H(r) the last two perturbations satisfy 



2^2 



r(r 4 — r^)(l — t> 2 ) r 5( r 4 _ r^)(l — v 

z ,, r) Zvrj _ 2vri (3r 4 (l - v 2 ) + r 4 (2*, 2 - 3)) 



(155) 



r 4 - 3r 4 (l - v 2 ) v V(r 4 - r 4 )(l - t; 2 )(r 4 - 3r 4 (l - v 2 )) ' 

H> (r) = - H(r) + ? (r MV 2 'j-2r*(l -v*)) 

) K ) r(r 4 -ri)(l-v 2 ) K } {rh\ - 3r 4 (l - i; 2 )) 2 

Hi \ ( ^ u 2r * (v 2 r A h il + v 2 ) + r\l + 2v 2 -^)) \ 

^^ r 4_ r 4+ r 5( r 4_ r 4)( 1 _ t ,2)2 ) + 

2r 4 (-9r 8 (l - t; 2 ) 2 (l + v 2 ) + u 2 rg(3i; 2 - 5) + r 4 r 4 (7 + v 2 - lh/ + 3v 6 )) 
+ W r(r 4 - r 4 )(l - i; 2 )(r 4 - 3r 4 (l - w 2 )) 2 + 

+ g V(l-, 2 ) 



(r 4_ r 4 )(r 4_ 3r 4 (1 _ ?; 2 )) • 

(156) 

We can now find the asymptotic behavior for all perturbations close to the horizon and close 
to the boundary. The results are summarized in table (A). 

The behavior of H 0i and Hu close to the boundary are consistent with the equations of 
motion for the boundary stress energy tensor. The last can be rewritten as 

5T m = , 5T n = , (157) 

where 5T^ U denote the perturbations away from the ideal boosted fluid stress energy tensor. 
From the AdS/CFT dictionary we know that 



5T n = lim r 2 # Ui (r) = , 5T W = lim r z H ll (r) = . (158) 

r— >oo r— >oo 

in agreement with (157). 



B Expansion near the boundary 

Given a conserved boundary energy momentum tensor the boundary conditions for an asymp- 
totic AdS metric are fixed. There is a procedure [67] that allows one to find such metric 
expanded in powers of - where r is the radial coordinate in the Poincare AdS patch, such 
that r = is the horizon and r = oo the boundary. In our case the procedure simplifies. We 
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r — )• oo 


r -+r h 


r — >■ 77 


Z{r) 




(r — 4r hV(i-v' 2 ) 


Co + Ci(r - rt) + C2(r - r/) 2 + • • • 


H 00 (r) 




qv i — 


/ l (°) + /l (°)( r _ r/ ) + ... 


H(r) 




y^ 2 {r-r h ) «r h v(i-^) jV < 

L V ' V 


/i + M r - rf) H 


n \r) 


-9 

T 


2ii)(-igD+2(l-2i> 2 )v / l--« 2 -(J 2 (l-f 2 ) 3/2 ) 


^0 + "1 V — r f ) ~r 


(3^ 2 -2)(^-4jVl-f 2 ) 

..• 91' i— 

■(r-r h ) 4r hV(i-« 2 ) ; i; < W| 


H i(r) 






C ) + M U1, (r-r / ) + --- 








/i 11) + ^ 11) (r-r / ) + --- 



Table 1: The behavior of the perturbations close to the boundary, the horizon and the 
pole of the equation (150) is given. Here, r/ = ^374 ^^1/4 is the location of the pole of 
(150). The boundary conditions for the perturbations at the boundary r — > 00 are that the 
metric is unchanged from the Minkowski metric. At the horizon the condition is that the 
asymptotic behavior corresponds to an infalling graviton in the unboosted black hole. In the 
linearized approximation the overall scaling factor does not appear and we only show the 
power law behavior. The normalization of the perturbations close to the horizon are relative 
to the normalization of Z(r). The normalization of Zir) is taken to be 1 for the factor 

multiplying (r — 77J 4r ^ (I "" 2) . For v > J | the relative coefficients C ii2 are not computed. 

The asymptotic behavior of ifoi? #11 is not shown for brevity but can be easily inferred from 
(155)-(156). 



fix the energy momentum tensor to be 

T tt = e{x), T tx = Ci, T XX = C 2 , 



1 

-e 
2 



(159) 



which is obviously conserved (d fl T IJ ' u = 0). The energy density e(x) has to be computed from 
the hydrodynamic equations or a guess can be made. In any case this fixes the boundary 
condition and allows us to extend the metric as: 
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where the boxed terms are fixed by the boundary conditions and the rest can be computed 
from solving the Einstein equations. In the absence of an exact metric for the shock wave, 
this expansion provides more information about it and could possible be used in the future as 
a check of given solutions and as a starting point for a numerical method. Although we show 
a few terms, it should be noted that using a computer algebra program we found easily the 
expansion up to order -§o although it is too lengthy to display here. These are enough terms 
to attempt a reconstruction of the metric using Pade approximants. The condition that 
determines the function e(x) then comes from demanding that the metric does not develop 
a singularity. In fact, in Fefferman-Graham coordinates the AdS-Schwarzschild black hole 
metric becomes degenerate at the horizon and one cannot go beyond the horizon in these 
coordinates. 
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